Theory of Spin Relaxation in Two-Electron Lateral Coupled Si/SiGe Quantum Dots 
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Highly accurate numerical results of phonon-induced two-electron spin relaxation in silicon double 
quantum dots are presented. The relaxation, enabled by spin-orbit coupling and the nuclei of 29 Si 
(natural or purified abundance), are investigated for all relevant parameter regimes, the interdot 
coupling, the magnetic field magnitude and orientation, and the detuning. We calculate all relaxation 
rates for zero and finite temperatures (fOO mK), concluding that all findings for zero temperature 
qualitatively remain valid also for fOO mK. We confirm the same anisotropic switch of the axis of 
prolonged spin lifetime with varying detuning as recently predicted in GaAs. However, there is a 
striking difference compared to the GaAs counterpart. In silicon, the hyperfine-induced relaxation 
rate is negligible in all cases we studied — even for natural silicon. The spin-orbit coupling, although 
weak, is the dominant contribution, yielding anisotropic relaxation rates of at least two order of 
magnitude lower than in GaAs. 

PACS numbers: 72.25.Rb, 03.67.Lx, 71.70.Ej, 73. 21. La 
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I. INTRODUCTION 

Since the proposal of Loss and DiVincenzo, 1 electron 
spins in semiconductor quantum dots have been in the 
perpetual focus of research on spintronics. 2,3 In GaAs 
based qubits, which are the state of the art, the essen- 
tial gate operations 4-7 for quantum computation 8,9 have 
been demonstrated. 10-14 But GaAs possesses a serious 
handicap for coherent spin manipulations — the nuclear 
spins. 15,16 Controlling this source of decoherence is of 
major interest and an active field of research. 17-22 

An alternative to III-V semiconductors with inherent 
nuclear spins are systems composed of atoms without 
nuclear magnetic moment, such as Si and C. 23-25 Natu- 
ral silicon consists of three isotopes: 28 Si (92.2%), 29 Si 
(4.7%), and 30 Si (3.1%). 26 Hereof only 29 Si has non- 
zero nuclear spin (1=1/2), and purification can further 
reduce its abundance down to 0.05%. 27 ' 28 For this rea- 
son, silicon-based quantum dots have become the new 
focus of interest, and recent progress emphasizes their 
perspectives. 29-32 Another advantage of silicon 9,25 over 
GaAs is a larger g factor, which allows spin manipula- 
tions in smaller magnetic fields. On the other hand, de- 
vice fabrication of silicon dots is more challenging, 33 the 
spin-orbit interactions are weaker, and the dots must be 
smaller due to a larger effective mass. 

Bulk silicon has six equivalent conduction band min- 
ima located on the A-lines, at k v 0.84fco toward the 
six X points of the Brillouin zone, where kg = 2ir/ao 
with ao = 5.4 A the lattice constant. 34-36 They are typ- 
ically referred to as A- valleys or X- valleys. In general, 
their degeneracy is lifted by strain, or by the presence 
of an interface. 35,37 In a (OOl)-grown silicon heterostruc- 
ture, the four in-plane valleys are split by at least 10 
meV from the two lower-lying ±z valleys, resulting in 
a twofold conduction band minimum. This remaining 
degeneracy is further split if the perpendicular confine- 
ment is asymmetric, resulting in an energy difference 



called the ground state gap. 38-45 As the valley degen- 
eracy is believed to be the main obstacle for silicon- 
based quantum computation 42,46,47 , a large valley split- 
ting is desired. If this is the case, the multi- valley sys- 
tem can be reduced to an effective single valley qubit, 
a potentially nuclear-spin-free analog to the well-know 
GaAs counterpart. 46,47 In fact, many recent experiments 
performed on Si/SiGe quantum dots have no evidence 
of valley degeneracy, 29-31,48-50 indicating that the split- 
ting is large enough to justify a single valley treatment. 
On the other hand, a recent proposal of valley-defined 
qubits uses the valley degree of freedom as a tool for gate 
operations. 51 This requires precise control of the ground- 
state gap, a challenging task for the future. In this work 
we assume that the valley splittings are larger than the 
typical energy scale of interest so that the effective single 
valley approximation is valid. 

The spin relaxation and decoherence has been investi- 
gated theoretically and experimentally in silicon-based 
single and double dots from single to many electron 
occupancy. 30,31,46,49,50,52-63 Our work completes these 
findings by a global, quantitative understanding of two- 
electron lateral silicon double quantum dots. We investi- 
gate the spin-orbit and hypcrfinc induced relaxation rate 
as a function of interdot coupling, detuning, and mag- 
nitude and orientation of the external magnetic field for 
zero and finite temperatures, and for natural and isotopi- 
cally purified silicon. We pay special attention to the spin 
hot spots 64 , and investigate individual relaxation chan- 
nels. We assume the double dot orientation with respect 
to the crystallographic axes that is used most often in 
experiments. This work is an extension to the findings 
in Ref. 65 for GaAs, and we highlight the differences be- 
tween these two materials. 

We find that due to the small spin-orbit coupling the 
spin relaxation rates are typically at least two orders of 
magnitude lower than in comparable GaAs dots, and that 
the relaxation rate peaks at spin hot spots are very nar- 
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row in parameter space. For detuned double dots, the 
energy spectrum close to the singlet-singlet anticrossing 
is qualitatively different from the GaAs counterpart, due 
to the rather small single-dot exchange coupling com- 
pared to the anticrossing energy. We also find that the 
hyperfinc-induccd relaxation rates of natural silicon are 
at least two orders of magnitude lower than the spin- 
orbit induced relaxation rates, even for the unpolarized 
triplet of detuned double dots contrariwise to GaAs. 65 
The hypcrfine- induced rates of purified silicon are, again, 
about two orders of magnitude lower compared to natu- 
ral silicon. We conclude that, concerning the relaxation, 
the nuclear field is negligible. Thus, the anisotropy of 
the spin-orbit field manifests in all relaxation rates we 
calculated, yielding the electrically controlled directional 
switch of the easy passage 66 (a particular orientation of 
the magnetic field for which the relaxation as a function 
of some parameter is significantly lower than for other ori- 
entations), previously found in GaAs. 65 A temperature of 
O.IK does not change our findings in any qualitative way. 
We give quantitative results for experimentally realistic 
parameters covering the whole experimentally significant 
range of the parametric space. 



II. MODEL 

We consider a z = [001] grown top-gated Si/SiGe het- 
erostructure defining a laterally coupled double quantum 
dot within the silicon layer with a fraction of 29 Si iso- 
topes. The double dot is charged with two electrons and 
not coupled to leads. Assuming the validity of the ef- 
fective single valley approximation, 58 the Hamiltonian in 
the two-dimensional and the envelope function approxi- 
mation reads 



H 



i=l,2 



(Tj + V t + H z>i + H SOtl + tf nuc>i ) + H c . (1) 



The operators of position r and momentum P are two- 
dimensional, where x = [100] and y = [010]. The single- 
electron terms are labeled by the electron index i. The ki- 
netic energy is T — P 2 /2m, with the kinetic momentum 
P = — iftV + eA, the effective electron mass m, and the 
electron charge — e. For an external magnetic field, given 
by B = cos 7, -By sin7,S z ), where 7 is the angle be- 
tween the in- plane component of B and x , the vector 
potential in symmetric gauge reads A = B z (— y,x) /2. 
We neglect the orbital effects of the in-plane magnetic 
field component, which is a good approximation up to 
roughly 10 T for common heterostructures. 67 The elec- 
trostatic potential, 
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2ml% 



min{(r 



d) 2 ,(r + d) 2 } + e E-r, (2) 



consists of the biquadratic confinement 47 ' 68 and the ex- 
ternal electric field. For E — 0, the potential is mini- 
mal at ±d, calling 2<i/Z the interdot distance. The sin- 
gle dot scale is given by the confinement length l , and 



equivalently by the confinement energy Eq = h 2 /(to/q). 
The electric field E is applied along the dot main axis 
d, where the angle 8 gives the in-plane orientation with 
respect to x. Turning on E shifts the potential minima 
relative to each other by the detuning energy e = 2eEd. 
The geometry is plotted in Fig. 1 of Rcf. 69. 

The Zeeman term is Hz = (5/2) ^bc-B, with the vec- 
tor of Pauli matrices er = (<r x , a y ,a z ), the effective Lande 
factor g, and the Bohr magneton p^- The spin-orbit cou- 
pling, H so = H br + iJ d , includes the Bychkov-Rashba 2 ' 70 
and the generalized Dresselhaus Hamiltonians, 2,71-73 

Hbr = (fr/2ml hl ) (<T x P y - <r y P x ) , (3) 
H d = (h/2ml d ) (-a x P x + a y P v ) , (4) 

parameterized by the spin-orbit lengths Ibr and l d . In 
this work we assume interface inversion asymmetry and 
choose Ibr and Id of comparable strength, according to 
Rcf. 73. The nuclear spins of 29 Si dominantly couple 
through the Fermi contact interaction 16,74 ' 75 



H nuc — p' S )^I n ■ a- <5(R — R n ), 



(5) 



where (3 is a constant, I„ is the spin of the n-th nucleus at 
the position R„, and R = (r, z) is the three-dimensional 
electron position operator. Here we need to consider the 
finite extension of the wavefunction perpendicular to the 
heterostructure interface. We assume it is fixed to the 
ground state of a hard- wall confinement of width w. This 
defines the effective width, 



h z = 



I 



dz\?P(z)\ 4 



(6) 



which evaluates to h z = 2u>/3. Finally, the Coulomb 
interaction is Hq = e 2 /Aire \ y\ — y?\, with the material 
dielectric constant e. 

The energy relaxation is enabled by phonons, whereas 
spin-orbit interactions allow for a spin-flip. In a (001)- 
grown quantum well of silicon, the electron-phonon cou- 
pling for intravalley scattering is the deformation poten- 
tial of transverse acoustic (TA) and longitudinal acoustic 
(LA) phonons, given by 56,60,76 - 79 



I HQ 
2pVc x 



^,Ae iQR -&Q,Ae- iQR l, (7) 



where 



£>£ = (S d e^ • Q + E u e^ z Q z ) 



(8) 



The phonon wave vector is Q = (q, Q z ), and Q = Q/Q. 
The polarizations are given by 80 A = TAl, TA2, LA, the 
polarization unit vector reads e, and the phonon annihi- 
lation (creation) operator is denoted by b (tf). The mass 
density, the volume of the crystal, and the sound veloci- 
ties are given by p, V, and C\, respectively. The efficiency 
of the electron-phonon coupling is set by the dilatation 
and shear potential constants, Sd and S u respectively. 
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tunneling energy T [meV] 

FIG. 1. (Color online) Calculated conversion from single- 
electron tunneling energy T to two-electron exchange coupling 
J neglecting nuclear spins. The scale for the interdot distance 
2d/lo is also shown. The cLITOW glVGS i/nuc) Eq. (10), of natural 
silicon. 

We define the relaxation rate (the inverse of the life- 
time T\) as the sum of the individual transition rates to 
all lower-lying states. Each rate (from \i) to is eval- 
uated using Fermi's Golden Rule at zero temperature, 3 

r « = £ t I^qI 2 l M ^l 2 - e q)> ( 9 ) 

where My = (i\e lC * -( R i+ R 2) |j) j s the matrix element of 
the states with energy difference Eij, and Eq is the en- 
ergy of a phonon with wave vector Q and polarization A. 
In this work we focus on the singlet (S) and the three 
triplets (T_,Tq,T + ) at the bottom of the energy spec- 
trum. 

Our numerical method is discussed in Refs. 81 and 
82. The extension to include the hyperfine coupling, 
Eq. (5), was introduced in Ref. 65. In this work, the two- 
electron basis for the configuration interaction method 
consists of 1156 Slater determinants, generated by 34 
single-electron orbital states. The discretization grid is 
typically 135 x 135. The relative error for energies is be- 
low 10 -5 . The reliability of our code is confirmed by the 
evaluation of Eq. (9) in an analytically solvable regime — 
weakly coupled dots in low magnetic fields. For details 
on this calculation, see the Appendix. 

We use the parameters of a SiGe/Si/SiGe quantum 
well grown along the z = [001] direction with a germa- 
nium concentration of 25%. The two-dimensional elec- 
tron gas is defined in the thin silicon layer with tensile 
strain. 34 ' 83 The in-plane effective mass is isotropic, given 
by the transverse mass of the X valley states, 37 and we 
use m — 0.198m e , 84 where m e is the free electron mass. 
The effective Lande factor is g = 2. 40 ' 85 Other mate- 
rial parameters read q = 9150 m/s (for LA phonons), 
Q = 5000 m/s (for TA phonons), p = 2330 kg/m 3 , and 
e = 11.9eo. 86 The choice of deformation potential con- 
stants is not unique, 78:79,87 and we use Sd = 5 eV and 
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FIG. 2. (Color online) Calculated energies of the lowest states 
varying (a) the interdot coupling (at B = 3 T), and (b) the 
detuning (at B — 0.5 T). Singlet states are given by dashed, 
triplets by solid lines. In (a), the energy of To is subtracted, 
and in (b), the quadratic trend in E is subtracted. The arrow 
in a) marks where J = E nuc (for natural silicon). 
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FIG. 3. (Color online) Schematic energy spectrum of an un- 
biased double dot showing the singlet S (dashed line) and the 
triplet To (solid line). For large interdot distances, the ex- 
change coupling J is given by the hyperfine splitting i?nuc, 
Eq. (10), and the eigenstates change to \]/[) and |ti). 



S u = 9 eV according to Ref. 86. The hyperfine coupling 
parameter reads (3 = — 0.1/xeV nm 3 , and 29 Si has spin 
7 = 1/2. For natural silicon, the 29 Si abundance is 4.7%, 
and we use an abundance of 0.05% for purified silicon. 
For the spin-orbit coupling strength we choose Zbr = 38.5 
/zm and id = 12.8 //m. 73,85 The confinement length is 
lo = 20 nm (Eq = 1.0 meV), in line with realistic dot 
sizes. 49,88 The double dot is oriented as d || [110]. The 
magnetic field is in-plane unless stated otherwise. 



III. RESULTS: UNBIASED DOUBLE DOT 

We parameterize in our model the coupling of the dou- 
ble quantum dot by the interdot distance 2d/lo- The 
corresponding experimental observables are the tunnel- 
ing energy T (single-electron occupancy) 58,89 and the ex- 
change coupling J (two-electron occupancy). 81 The con- 
version between these three equivalent parameters is plot- 
ted in Fig. 1 for clarity. 

The numerically calculated energy spectrum of the un- 
biased double dot is shown in Fig. 2. In this section, we 
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FIG. 4. (Color online) Calculated relaxation rates of (a) the 
singlet, (b) the triplet To, and (c) the triplet 7+ as a function 
of the in-plane magnetic field orientation 7 (angle) and the 
interdot distance 2d/lo (radius of the polar plot), for a double 
dot at B = 3 T. The x and y axes correspond to crystallo- 
graphic axes [100] and [010], respectively. The dot orientation 
d || [110] is marked by a line. The x axis is converted to the 
tunneling energy T and the exchange J, in addition to 2d/£o- 
The rate is given in inverse seconds by the color scale. The 
system obeys Ci v symmetry, so point reflection would com- 
plete the graphs. 

choose a magnetic field of B = 3 T. For the single dot 
(d = 0) the exchange coupling, J = E(T ) — E(S), is 
Jsd = 0.19 mcV. The Zeeman energy, Ez — g^B, ex- 
ceeds J for magnetic fields beyond 1.7 T. Consequently, 
we find in Fig. 2a that T_ is the ground state for all inter- 
dot distances. The singlet therefore has an anticrossing 
with an excited triplet in the strong coupling regime, here 
at J = 75 fieV for our choice of parameters. This scenario 
is hardly met in comparable GaAs double quantum dots, 
because the required magnitude of the magnetic field is 
above 10 T. 65 The silicon spectrum resembles the GaAs 
spectrum for magnetic fields below 1.7 T. 

At large interdot distances, the hyperfine coupling in- 
duces a splitting of S and To, given by 



E nuc = 2 \ (^aTo|# nU c,i|<£ s S)| 2 . (10) 

In this regime, the lowest eigenstates are |f|) — (S + 
T )/J2 and || T ) = (S - T )/^2 (c.f. Fig. 3). We eval- 
uate Eq. (10) by averaging over random nuclear spin 
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FIG. 5. (Color online) Calculated channel resolved relaxation 
rates vs. interdot distance for both parallel (top) and per- 
pendicular to d (bottom) in-plane magnetic field orientation 
(B — 3 T). The relaxation channels of T+ are black, of To are 
blue, and of S are red. 



ensembles, and obtain E nuc w 2 neV for natural, and 
E nuc w 0.1 neV for purified silicon. This implies a 
crossover to the nuclear dominated regime at 2d/lo > 4.6 
(red arrows in Figs. 1, 2a, and 3) for natural, and at 
2d/lo > 5.3 for purified silicon. 

We plot the relaxation rates of the states S, To, and 
T + , denoted in Fig. 2, as a function of the interdot dis- 
tance and in-plane magnetic field orientation in Fig. 4. 
We also give the relaxation rates of individual channels 
for the two principal axes, that is for the in-plane mag- 
netic field component parallel and perpendicular to the 
dot main axis d, in the upper and lower panel of Fig. 5, 
respectively. We find that the relaxation rate of the 
singlet is highly anisotropic due to the effective, spin- 
orbit induced magnetic field. 65 ' 66 The rates are minimal 
if B || d, reaching the order of tens of milliseconds for any 
dot coupling strength (Fig. 5). We call this characteris- 
tic an easy passage. 66 In the strong coupling regime, the 
rate away from the easy passage is enhanced by orders of 
magnitude. This results from the coupling of the singlet 
with the excited triplet, which favors the transition into 
T_. For B || d, the rate at the anticrossing is extremely 
sensitive to variations of 7, such that the easy passage 
becomes very narrow. A spectral crossing, which is lifted 
in the presence of spin-orbit coupling, is denoted a spin 
hot spot. 64 

The relaxation rate of To is given in Fig. 4b. We find 
the same general anisotropic behavior, which is that the 
rate is minimal for B || d. Figure 5 shows that the domi- 
nant channel of the relaxation is the transition To — ^ T_ . 
Consequently, there is no impact from the singlet-triplet 
anticrossing. However, the anticrossing of the excited 
triplet with Tq itself manifests in a very sharp peak of 
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magnetic field [T] 

FIG. 6. (Color online) Calculated relaxation rates of (a) the 
singlet, (b) the triplet To, and (c) the triplet T+ as a func- 
tion of the in-plane magnetic field orientation 7 (angle) and 
the magnetic field magnitude (radius of the polar plot), for 
a double dot with T = 0.1 meV. The layout with respect to 
the crystallographic axes is the same as in Fig. 4. The rate is 
given in inverse seconds by the color scale. 



its rate. This spike is also anisotropic, with a difference 
of roughly one order of magnitude (not visible in Fig. 4b 
due to its resolution). 

Panel c) of Fig. 4 shows the relaxation rate of T + . In 
addition to the anisotropic background, there are two 
spikes of enhanced rate generated by the anticrossings of 
T+ with the excited triplets. The enhancement close to 
the single dot regime originates from the dominant T + — > 
S transition (c.f. Fig. 5). Interestingly, the anticrossing 
of the singlet hardly influences the overall trend of this 
relaxation channel. 

We plot in Fig. 6 the relaxation rates of a weakly cou- 
pled double dot as a function of in-plane magnetic field. 
Here we find the same qualitative behavior for all three 
panels. Similarly as in Fig. 4, the relaxation rate is min- 
imal for B |j d, but there are no spin hot-spots here. 



IV. RESULTS: BIASED DOUBLE DOT 

In this section we consider a weakly coupled double 
dot with a finite detuning energy e. Figure 7 introduces 
important characteristic energies in a schematic energy 
spectrum. The state charge character is given in brack- 
ets: (1,1) indicates that there is one electron in each dot, 




FIG. 7. (Color online) Schematic energy spectrum of a biased 
double dot without magnetic field. The singlets are given in 
dashed lines, the triplets in solid lines. 



and (0,2) states that both electrons are in the same dot. 
In the spectra, we subtract the quadratic trend in the 
electric field E. This way, the (1,1) states are displayed 
horizontally unless influenced by anticrossings. The im- 
portant quantities are the single-dot exchange coupling 
Jsd, the double-dot exchange coupling J, and the sin- 
glet and triplet anticrossing energy splittings, labeled in 
Fig. 7 respectively. The single-dot exchange coupling Jsd 
is set by material parameters, i.e. the Coulomb interac- 
tion, and system parameters, i.e. the confinement length. 
For non-interacting electrons, Jsd is equal to the confine- 
ment energy Eq, here 1 meV. For interacting particles, 
the Coulomb repulsion has strong impact on the symmet- 
ric ground state, the singlet, as here the electrons tend 
to group together. The first excited state, the triplet, is 
antisymmetric with respect to point reflection at the dot 
origin, and therefore less affected. As a consequence, Jsd 
decreases as the Coulomb interaction strength increases. 
For our choice of parameters Jsd = 0.2 meV. In con- 
trast, Jsd increases as the confinement length decreases. 
For instance, a confinement length of Iq = 17 nm results 
in Jsd ~ 0.3 meV. This can be understood as follows. 
On the one hand, a stronger confinement increases the 
Coulomb strength due to smaller effective particle dis- 
tances |ri — r2 1 in Hq. This is an effect somewhat linear 
in Iq . Then, one could expect Jsd to decrease. How- 
ever, the confinement energy Eq scales as Iq 2 , by what 
the exchange coupling increases in a similar way. This 
scaling dominates, such that the single-dot exchange cou- 
pling increases. The double-dot exchange coupling J de- 
creases exponentially with increasing interdot distance 
(^ 58,89 j n weak coupling regime it holds J <C Jsd, 
and we choose d such that J = 0.6 fieV. The anticross- 
ing gap of a spin-alike pair of states at the (1,1) «-> (0,2) 
transition depends on the interdot distance as well. For 
increasing d (decreasing J), these gaps decrease, that is 
the anticrossings vanish as 2d/lo — > 00. 

The numerically calculated energy spectrum is plot- 
ted in Fig. 2b for a magnetic field of B = 0.5 T. The 
spectrum is qualitatively different from the GaAs double 
dot counterpart (see Fig. 1 in Ref. 65). In a compara- 
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FIG. 8. (Color online) Calculated relaxation rates of (a) the 
first excited state (S or T_, see Fig. 2b), (b) To, and (c) T+ as 
a function of the in-plane magnetic field orientation 7 (angle) 
and detuning energy (radius of the polar plot), for a double 
dot with 2d/l = 2.85 (T = 0.1 meV), and B = 0.5 T. The 
layout with respect to the crystallographic axes is the same 
as in Figs. 4 and 6. The rate is given in inverse seconds by 
the color scale. 

ble GaAs double dot, the singlet and triplet anticrossings 
gaps are small compared to the single-dot exchange cou- 
pling. Consequently, the singlet anticrossing is well sep- 
arated from the triplet anticrossing, and the excited sin- 
glet is close to To in between these anticrossings. In this 
regime, the relaxation rate is dominated by the coupling 
between T and the excited singlet, and this coupling is 
solely due to the nuclear spins. 65 In silicon (c.f. Fig. 2b) 
the anticrossings are much more pronounced, and Jsd is 
smaller than in GaAs. Therefore, in the experimentally 
most relevant detuning energy range up to the triplet- 
triplet (anti-)crossings, the excited singlet is far away in 
energy. Thus, there is no significant (hyperfine) coupling 
between To and the excited singlet, i.e. the relaxation 
rates are not affected by the excited singlet in any case 
(see below). 

We plot the relaxation rates of the detuned double dot 
in Fig. 8. Panel a) corresponds to the first excited state, 
that is S for detuning energies up to 1.97 meV, and T_ 
beyond. At the singlet-triplet anticrossing, the relax- 
ation rate is very low as the transferred energy becomes 
very small. The effective, spin-orbit induced, magnetic 
field, 65,66 Eq. (12), confirms that the easy passage occurs 
if the external, in-plane magnetic field is perpendicular 
to d. The anisotropy is also visible for the relaxation 
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FIG. 9. (Color online) Calculated relaxation rates of a biased 
double dot in an in-plane magnetic field (B — 0.5 T, 7 = 
37r/4) as a function of detuning. The straight lines give the 
spin-orbit induced relaxation, the wiggly lines the hypcrfinc 
induced relaxation rates (natural silicon). 

rates of To and T + , panel b) and c) respectively. This di- 
rectional switch has also been found numerically in GaAs 
dots. 65 There is no signature of the singlet-triplet anti- 
crossing in the rate because of the exact compensation of 
individual relaxation channels. 65 Also, there is no indi- 
cation of the crossing of T + with the excited triplet T_ . 
Other anticrossings with excited triplets (at e rj 2.47 
meV) manifest in extremely narrow peaks of the rate, 
not visible in Fig. 8 at the current resolution. 

Let us comment on the role of hypcrfinc spins. Com- 
paring the interaction strengths with the spin-orbit fields, 
they are expected to be negligible. Indeed, the Over- 
hauser field characterizing the fluctuating collective nu- 
clear field, 16 




of natural silicon is of the order of tens of /iT, for pu- 
rified silicon even one order of magnitude lower. On 
the other hand, the effective spin-orbit field (xd, yd — 
(x±y)/V2), 66 

B so - Bx{x d (C -£ X )[1T0] +y d (C +Z7 1 )[110]}/V2, 

(12) 

is about 2 mT at B — 1 T. Still, in GaAs we have found 
that despite a similar discrepancy, there are cases where 
the nuclear field dominates the spin-orbit field, as the lat- 
ter is quenched by symmetry imposed selection rules. 65 
Here, however, such anomalous situation does not arise 
due to the different spectrum, as we discussed above in 
Sec. IV. The matrix elements of the spin flipping tran- 
sitions are in Si therefore always dominated by the the 
spin-orbit fields, rather than nuclear spin fields and the 
same holds for anticrossing gaps. An illustration is given 
in Fig. 9. 
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The second possibility we considered, was a dot de- 
tuned so far (so small J), that singlet and triplet T be- 
come degenerate with respect to ii^uc- The Hamiltonian 
eigenstates change, from entangled states into separable 
states with spin up or down in the left or right dot re- 
spectively (we show this schematically in Fig. 3). Figures 
in Sec. Ill cover this regime but the qualitative change 
in the eigenstate character has no visible effects on the 
relaxation rates (verified also for 2d/lo > 5, not shown). 

Next we considered direct transitions due to random 
nuclear fields without phonon assistance. Such transi- 
tions are possible if the eigenstates have unsharp energies 
(finite lifetimes). As the states we are interested in are 
low lying, even at finite temperature their energy broad- 
ening is so small that the resulting nuclear induced spin 
relaxation is negligible. 

Finally, we considered the consequences of the random 
character of the nuclear field, which blurs the electron en- 
ergies. This statistical, rather than quantum mechanical, 
uncertainty can be grasped roughly by convoluting the 
relaxation curves with a Gaussian with an appropriate 
width, depending on which parameter we change, defined 
ultimately by the energy g/iBB nuc . We find this width 
to be unnoticeable small — as an example, the extremely 
narrow peaks in Fig. 5 survive practically untouched by 
such smoothening. We conclude that unpolarized nuclear 
spins in natural or purified Si are not expected to be vis- 
ible in any way in the electron relaxation. 

The figures presented and results discussed in this arti- 
cle were for zero temperature. In our model, a finite tem- 
perature amounts solely to allowing for energy increasing 
transitions (phonon absorption), in addition to phonon 
emission processes only which are present at zero tem- 
perature. We analyzed this possibility, adopting a typ- 
ical experimental value of 100 mK. We have not found 
any case where such additional transitions would change 
the relaxation rates in any significant way (figures not 
shown) . Our conclusion from these investigations is that 
the relaxation character, most notably its anisotropics, 
will not be influenced by experimentally relevant sub- 
Kelvin temperatures. 
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Appendix A: Analytical Calculation of Relaxation 
Rates 

In this section we calculate the relaxation rate, Eq. (9), 
analytically adopting several approximations. The cal- 
culations prove useful to explain the physical mechanism 
and to verify our numerical results. The validity of the 



approximations will be discussed afterwards. In the fol- 
lowing, the hyperfine coupling is neglected. 

Approximating the sum in Eq. (9) by an integral, and 
rewriting the (^-function with respect to the z component 
of Q, we obtain (i ^ j) 

x [l(Q.-Qi) + S(Q,+Qi)] , (Al) 

where Q\ = E 2 -/(h 2 c\) — q 2 . Assuming the validity of 
the dipole approximation, the matrix element reads 

Aftf«i<i|q-(n+r 2 )|.7>, (A2) 

where |i) and \j) are the spin-orbit coupled two elec- 
tron eigenstates. Note that the contribution of the wave 
function overlap along the z direction in My is about 
l, 90 which is consistent with the two-dimensional ap- 
proximation. We restrict ourselves to weakly coupled 
double dots, i.e. d ^> Iq, and incorporate the effect of 
spin-orbit coupling perturbatively via a Schrieffer- Wolff 
transformation. 91-93 The eigenstates then read (I labels 
the electrons) 

w-- i0 fi^EEf^f4 (A3) 

with the transformation operator O = — J2i n so,l • o"j/2, 
where 

( Xl Vl Xl yi n\ (\a\ 
n so,i = I 7 -j—, -j -r, , (A4) 

\<d V V I'd J 

and the effective spin-orbit operator H so ,i = H^ Q i+H^\, 
where 

H^,i = f /«B(n scM x B) • a u (A5) 

^ = ±^-^L s ^+ const. (A6) 

Here, L z = l z + (e/2)r 2 B z , where l z is the operator of 
angular momentum. The states in Eq. (A3) labeled with 
subscript are eigenstates of the Hamiltonian 

H° = ( T i + V i + H z,i) + H c , (A7) 

i=l,2 

their eigenenergies are denoted as E°. We use the 
Heitler-London ansatz 47 to approximate the eigenstates 
of Eq. (A7). 

We use Eq. (A3) to evaluate the matrix element My. 
It is straightforward to show that contributions from 
coupling within the lowest four-dimensional subspace, 
M. = {S,T-,Tq,T + }, are zero or exponentially sup- 
pressed in d/lo. As a result, the relaxation requires cou- 
pling via higher states. Neglecting the L z contribution 
to the effective spin-orbit coupling, Eq. (A6), wc obtain 
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Ma 



EE 

k^M 1 = 1,2 



o(i|(n SOj ; x B) • cri\k) 



a{k\q x xi + q y yi\j)n 



)<fc|(n SOj ; x B) • (Ti\j} 



o{i\q x Xi + q y yi\k) 



(A8) 



The singlet is symmetric with respect to the inversion 
operator / (point reflection in real space) , the triplets are 
antisymmetric. 81 Consequently, it follows from Eq. (A8) 
that, within the dipole approximation, the singlet-triplet 
transition is forbidden. We also find that Eq. (A8) forbids 
a T + -H> T_ transition because the effective spin-orbit 
operator, H^ a ; , acts on only one of the two electron spins. 
Let us now look at the transition between To and T± . 

To evaluate Eq. (A8), we reduce the infinite sum over 
k to cover only states within the energy window of about 
the confinement energy, h 2 /(ml^). Additionally, we can 
exclude any singlet from the sum, because the electron- 
phonon operator does not act in spin space. What is left 
can be captured by the Heitler-London approach. 

Let |R0) be the (orbital) ground state of a single dot 
shifted to the "right" by d, i.e. the Fock-Darwin state of 
the right dot with the principal quantum number n = 
and the orbital quantum number I — 0. Analogously we 
define the ground state of the "left" dot. The properly 
symmetrized triplet lowest in energy is 



(A9) 




|tf T ) = (|R0,L0) - |L0,R0)) <g> |T)/V2 



in-plane magnetic field [T] 



FIG. 10. (Color online) Calculated relaxation rates of indi- 
vidual transition channels as a function of in-plane magnetic 
field for a weakly coupled quantum dot. The magnetic field 
is orientated along [110] (7 = 37r/4), and the dots along [110] 
(8 = 7r/4). The interdot distance is 2d/lo — 2.85, yielding the 
tunneling energy T = 0.1 meV. The dashed, black line gives 
the analytical relaxation rate, evaluated with Eq. (A14). 



The orbitally excited triplets can be constructed analo- 
gously, using |R1), and |L1), the displaced Fock-Darwin 
states with n — and |^| = 1: 

\k±)o =(|R0,L1) - |L1,R0) ± (|R1,L0) - |L0,R1))) 

<g> |T)/2. (A10) 

Neglecting the wave function overlap of states local- 
ized in different quantum dots, we calculate the ma- 
trix elements o(^t|^i \k±)o and o(^T\yi\k±)o analyti- 
cally, yielding 



{9 T \xi\k^ 



I0/V8, 

sgn(Z)y /v / 8, 



(All) 
(A12) 



and {V T \xi\k-)o = o^tMMo = 0. We use 
Eqs. (All) and (A12) as an approximation for the ma- 
trix elements in Eq. (A8). We also require the matrix 
elements of Pauli matrices respecting the spin quantiza- 
tion axis along B. They read 



2V2 



'cos (7 
sin(7 



(T ± \tr l \T ) = 

- 9) + cos(7 + 1 

- 0) + sin( 7 + ( 

2 sin(0) 



) ± 2isin(7) N 
=F 2i cos(7) 



(A13) 



where 9 = cos 1 (B Z /B^). The energy differences in 
Eq. (A8) are approximated by the confinement energy, 
tf/{mll). 



With these ingredients, we can solve Eq. (Al), inte- 
grating over the phonon momentum, and obtain 



- f„-7 ( 3 cj2 1 



12-nph 10 

-~1 I 3 ~2 1 2~ 77 , ~2\ _ 
A \yi V35"u ~t~ 5" u " d + "dl ' "t 35"u 

with the effective spin-orbit length C so , defined by 



C ,- 7 4S?,1 , (A14) 



£ 2 = 



2 (C 



r d 2 ) 



sin(27) 



if 9 = 0, 
if 9 = tt/2. 



(A15) 



Now we discuss the validity of the approximations 
used during the derivation of Eq. (A14). The matrix 
element Mjj is calculated using the dipole approxima- 
tion, Eq. (A2). It requires that the energy difference 
between the transition states, here To and T±, fulfills 
Eij -C Hcx/Iq. 90 Using J3y = g^B, and c\ — ct, we ob- 
tain the condition B <C 1.4 T. We consider also a weakly 
coupled double dot, d 3> lo, to comply with most ex- 
periments. This limit ensures negligible matrix elements 
among the states of Ai, and justifies the Heitler-London 
approximation. Here, the spectrum also develops bun- 
dles of eigenenergies separated by the confinement energy 
H 2 /(mlQ), a fact used to approximate the energy differ- 
ences in Eq. (A8). Note that within the restriction of the 
dipole approximation (B <C 1.4 T), the Zeeman energy 
(Ez <C 0.16 meV) is negligible compared to the confine- 
ment energy (Eq = 1 meV). The Schrieffer- Wolff trans- 
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formation is the essential tool for a perturbative treat- 
ment of spin-orbit coupling in the double dot. 89 Perturba- 
tion theory with the unitarily transformed Hamiltonian 
yields results which are higher order in small quantities 
compared to the original Hamiltonian. Finally, we note 
that, since L z is symmetric with respect to the inver- 
sion /, the perturbation ; , Eq. (A6), vanishes for all 
transitions To -R- T±. 

We compare the analytical formula for the relaxation 
rate of the transition T — » T_, and T + — > T , given in 
Eq. (A14), with the numerical results in Fig. 10. We find 



perfect agreement for low magnetic fields, in line with the 
condition B <C 1.4 T. For larger magnetic fields, the re- 
sults significantly deviate deviate from the B 7 power law, 
due to the break down of the dipole approximation. We 
also find that the S -B- T± relaxation channels, which we 
found to be zero in the lowest order dipole approximation 
due to their symmetry, show B 9 dependence, indicating 
that the relaxation is driven by the second order term of 
q. Being of higher order, the relaxation rate, for small B, 
is at least one order of magnitude lower than the T o T± 
transitions. 
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